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ABSTRACT 

Aims. Bow shock waves are a common feature of groups and clusters of galaxies since they are generated as a result 
of supersonic motion of galaxies through the intergalactic medium. The goal of this work is to present an analytical 
solution technique for such astrophysical hypersonic blunt body problems. 

Methods. A method, developed by Schneider (1968, JFM, 31, 397) in the context of aeronautics, allows calculation of 
the galaxy's shape as long as the shape of the bow shock wave is known (so-called inverse method). In contrast to 
other analytical models, the solution is valid in the whole flow region (from the stagnation point up to the bow shock 
wings) and in particular takes into account velocity gradients along the streamlines. We compare our analytical results 
with two-dimensional hydrodynamical simulations carried out with an extended version of the VH-1 hydrocode which 
is based on the piecewise parabolic method with a Lagrangian remap. 

Results. It is shown that the applied method accurately predicts the galaxy's shape and the fluid variables in the post- 
shock flow, thus saving a tremendous amount of computing time for future interpretations of similar objects. We also 
find that the method can be applied to arbitrary angles between the direction of the incoming flow and the axis of 
symmetry of the body. We emphasize that it is general enough to be applied to other astrophysical bow shocks, such 
as those on stellar and galactic scales. 
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1. Introduction 

Astrophysical bow shocks are ubiquitous in the Universe 
and can be observed on all scales, from the Earth's bow 
shock through the heliospheric (driven by the moving solar 
wind) and on to corresponding stellar wind bow shocks. On 
even larger scales, we see bow shocks of galaxies in the in- 
tergalactic medium (IGM) in groups and clusters, and pos- 
sibly even in mergers of galaxy clusters ( |Markevitch et ah 
2002 ) . The main ingredients of the feature are a compress- 
ible medium and a body moving through it supersonically. 
This was already investigated decades ago in aeronautical 
engineering in the context of supersonic aircrafts, and is 
commonly known as the supersonic blunt body problem. 
However in astronomy, with the exception of a few papers, 
analyt ical solutions have bee n scarce. For example, the pa- 
per of Canto & Raga (1998]) considers the supersonic mo- 
tion of a spherical body and treats the post-shock flow in 
a thin-shell approximation, which unfortunately restricts 
the applicability of the method to more realistic problems. 
The purpose of this paper is to introduce a more general 
analytical method, which allows calculating the stand-off 
distance of the shock and the complete post-shock flow. To 
demonstrate this, the analytic solutions will be compared 
to numerical two-dimensional simulations. It will be shown 
that the method even works reasonably well in the case 
of low Mach number shocks. As an important example of 
great astrophysical relevance, we analyse the bow shock of 
a galaxy moving supersonically through the IGM and show 



that the associated soft X-ray emission can be modelled 
fairly well. 

About half of the galaxies in the Universe are found in 
groups and clusters, which are large complexes of galax- 
ies held together by the mutual gravitational attraction of 
their members, intergalactic gas, and, above all, dark mat- 
ter. With the advent of imaging telescopes in X-ray astron- 
omy, it became evident that galaxy clusters are intensive 
sources of X-ray radiation owing to the hot plasma located 
between the galaxies and trapped in the group's/cluster's 
potential well. Because of their low relative speeds, galax- 
ies in groups and poor clusters affect each other and the 
surrounding IGM gravitationally stronger than the faster- 
moving cluster galaxies, leading to a variety of fascinating 
interaction or even merging processes. Galaxies dashing at 
super- or hypersonic velocities through the IGM produce 
wakes of gravitationally focused gas, lose mass due to ram 
pressure, thereby injecting metals into their gaseous envi- 
ronment. Finally, leading bow shock waves are generated 
ahead of the galaxies that change the state of the gas ir- 
reversibly and thus play an importa nt role in the galaxi es' 
structure and subsequent evolution (Stevens et al. II 1999). 

If 

ofl 



Trinchieri et al. 
plex X-ray emission 



2003 2005) have studied the com- 



the compact galaxy group Stephan's 
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Quintet (SQ) quite recently using Chandra and XMM- 
Newton observations. The prominent shock situated in the 
galaxy system has been resolved with Chandra into a nar- 
row north-south, somewhat clumpy structure between NGC 
7318ab and NGC 7319, which is more sharply bounded on 
the west side (probably due to a contact surface stabilized 
by a magnetic field that is indicated by radio continuum 
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Fig. 1. Shock-oriented coordinate system of boundary- layer 
type, where z and r are Cartesian coordinates for plane 
flow, x and y are the distances along the shock surface and 
normal to it, with u and v as the corresponding velocity 
components, Uqq is the free- stream velocity, and ip is the 
stream function. Flow quantities immediately behind the 
shock in the point N are denoted by a h at ( ~ ) , and in th e 
point S by an asterisk (*) (adapted from Schneide"r||1968 ) . 



emission) and embedded in a more extended diffuse emis- 
sion, which presumably represents preexisting IGM heated 
up by previous collisions (presuma bly NGC 7320c was in- 



volve d) . In the simplest scenario, as |Trinchieri et aL (2003 
2005) report, the shock results from the high- velocity col- 



ision of the gas-rich spiral galaxy NGC 7318b with previ 
ously stripped Hi gas in SQ. For an upstream Hi tem- 
perature of 100 K (and number density of about 6.5 x 
10 -3 cm -3 ), the velocity of the intruding galaxy is highly 
hypersonic (~ 1400kms _1 ), resulting in an upstream Mach 
number in the galaxy's rest frame of ~ 930. The gas 
inside the bow shock is then heated up to an (observed) 
energy of 0.5 keV. To explain these low post-shock temper- 
ature, the authors suggest an oblique shock scenario with 
a shock inclination angle j3 of about 30°. Oblique shocks 
are a general feature of bow shocks, which are ultimately 
responsible for the observed IGM structures. The analytic 
calculation of general bow shocks, however, is fairly compli- 
cated. Previous attempts bear severe restrictions, like the 
spherical shape of th e body and the t hin-sh ell approxima- 
tion of the model by | Cant 6 fc Raga| (JT998). On the other 
hand, bow shocks are ubiquitous phenomena in the inter- 
galactic and interstellar medium (ISM) , and an analytic de- 
scription like the one given in the present paper is therefore 
most desirable. 

The paper is organized as follows. In Section [2] the an- 
alytic solution for the blunt body problem is presented. In 
Section[3j this method is applied to Stephan's Quintet, and 
the results are also compared to numerical two-dimensional 
simulations. Section|4]closes the paper with our conclusions. 



2. Solution for the hypersonic blunt body problem 

A body moving supersonically through a gaseous medium 
will generate a shock wave, which is generally detached and 
curved. It appears to be planar at the body's nose and be- 
comes progressively more oblique sideways, until it decays 
into a sonic wave at large distances. As a result, this single 
shock wave represents all possible oblique shock solutions 
for the given Mach number with the wave angle rang- 
ing from (3 = 7r/2 to f3 = a, where a = arcsin(l/M) is 
the so-called Mach angle, i.e. the angle forming the cone to 
which small disturbances in a supersonic flow are confined. 
The blunt body shock layer, the volume between the body 
and the shock, is a mixed subsonic-supersonic flow, where 
the subsonic and supersonic regions are divided by sonic 
lines, i.e the loci of points at which the downstream Mach 
number is unity. Behind the normal, and nearly normal por- 
tions of the shock wave, the flow is subsonic, whereas be- 
hind the more oblique portions of the shock wave the flow 
is supersonic. It is precisely this hybrid nature of the flow 
that makes the blunt body problem extremely challenging 
mathematically, since the governing highly nonlinear par- 
tial differential equations of hydrodynamics are of elliptic 
type in the subson ic region and o f hyperbolic type in the 
supersonic region (|Anderson|2006). 

by Sch 



The analysis developed by Schneider I (1968) allows 



very elegant treatment of the inviscid hypersonic blunt 
body problem. It is an inverse method, which means that 
the shock wave shape is assumed, and both the body shape, 
which supports the assumed shock, and the flow field be- 
tween shock and body (i.e. the shock layer) are calculated. 
The fundamental advantage of this method over the ones by 
other authors is its uniform validity in the whole flow field 
(from the stagnation region up to large distances from the 
projectile nose). Until now this method has been adopted 
in aerospace engineering, in particular for solving the reen- 
try problem of space probes, space shuttles, etc., in plane- 
tary or terrestrial atmosphere. However, to the best of our 
knowledge, this paper demonstrates its first application to 
an astrophysical problem. 

We consider the plane or axisymmetric flow around a 
body and introduce a shock-oriented curvilinear coordi- 
nate system of boundary-layer type, where x is the distance 
along the shock surface in the plane formed by the shock 
normal and the direction of the uniform fluid flow, and y 
is the distance normal to the shock surface (Fig. HI. The 
corresponding velocity components are denoted by u and 
and z and r are the Cartesian coordinates for plane flow 
or the cylindrical coordinates for axisymmetric flow. The 
z-axis may be parallel to the direction of the incident flow. 
We emphasize that this is only assumed for convenience, 
but is no general restriction (see Section [4]). It can be de- 
duced from Fig. [I] that 

z = z + y sin P , (1) 



r = r — y cos f3 , 



(2) 



where $ is the shock inclination angle in the point N(z, f), 
with the shock-normal through Q intersecting the shock 
surface. Moreover, S is the point where the streamline 
through Q crosses the shock wave. The flow quantities im- 
mediately behind the shock in the point N are denoted by a 
hat ("), and in the point S by an asterisk (*). Undisturbed 
flow quantities far upstream are denoted by the subscript 
oo. 
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Since the functions z(x) and r(x) are known for a given 
shock shape, equations ([TJ and ([2| may be used to calculate 
the coordinates z and r of a point Q from its coordinates x 
and y. The curvature of the shock contour in the point N 
is denoted by k(x), defined as positive when the surface is 
concave on the side of positive y (cf. Fig. [I]) . The curvature 
of any of the other surfaces of constant y is U~ 1 



u = 1 - ky > . 



(3) 



The metric for this coordinate system is (cf. Hayes fc| 
Probstein||1966 ) 



ds 2 = H 2 dx 2 + dy 2 



(4) 



Thus, the governing hydrodynamical steady-state equa- 
tions for mass, momentum, and energy become 



dri pu dUr^ pv 



dx 



du 
dx 
dv 
dx 
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KU 



as 

dx 



dy 
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p dx 
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dS 
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(5) 


o, 
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o, 


(7) 


o, 


(8) 



where p is the fluid density, P the fluid pressure, S the 
entropy, r(x,y) the distance from the axis. The parameter 
j is or 1 for plane or axisymmetric flow, respectively. We 
used the entropy equation in place of the energy equation. 
We now want to make a von Mises transformation such that 
the independent variables (x, y) are replaced by (x = x, ip), 
where ^ is a stream function defined in terms of its partial 
derivatives 



dip 
dx 
dip 

dy 



(1 — ky)r 3 pv , 



-r 3 pu . 



(9) 
(10) 



The stream function is constant along a streamline and rep- 
resents the mass flow between the streamline ip = const, 
and the surface of the projectile, per unit depth for plane 
flows, and per unit azimuthal angle (in radians) for axisym- 
metric flows; i.e., ip satisfies the continuity equation. In the 
point TV, the stream function thus is simply 



1+3 



(11) 



In contrast, ip is connected with the coordinate r* of the 
point S by 



^ PocUoc-^ 



1+3 



In the new coordinate system with the variable y elimi- 
nated, the equations of energy, entropy, and momentum 
conservation read as 



2ft 



2ft* = const. . 



dS_ 

dx 



= 0, or S = S*(tl>), 



(13) 



(14) 



du dv 1 dP 

dx dx p dx 



= 0, 



(1-ky) (l-j|cos/3) 



,dP A 
r J — — - = km - 
dip 



dv 
dx 



(15) 
(16) 



where ft is the specific enthalpy. 

The variable y is a dependent variable after the von 
Mises transformation, and it obeys the equations 



dy_ 
dx 
dy_ 

dip 



(1 - %)' 



1 



{1 — j(y/r) cos /3}ri pu 



(17) 
(18) 



At this point it is useful to introduce the Landau sym- 
bol 0, which describes an asymptotic upper bound for the 
magnitude of a function in terms of another, usually sim- 
pler function; e.g., f(x) = 0(g(x)) means that \f(x)\ is not 
very large in comparison with |#(x)|. The flow quantities 
immediately behind the shock may be obtained from the 
Rankine-Hugoniot jump conditions in terms of the inverse 
compression ratio across the shock x = Pool P- They are 



u = /Too cos $ •> 
v £/ooXsin/3, 

P = Poo + PocUl(l - x) sin 2 /3, 



fto 



1 



Ul(l-x 2 )sm 2 p. 



(19) 
(20) 
(21) 

(22) 



These equations maintain their validity if the hats are re- 
placed by asterisks. The pressure within the shock layer 
is given by the Newton- Bus emann pressure law (see e.g. 
Hayes & Probstein 1966), a simplified form of the momen- 
tum equation (Il6[), which reads as 



P = P 



f3 



udtp' . 



(23) 



The following method is based on two main assump- 
tions. First, it is assumed that the inverse compression ratio 
across the shock is very small; i.e., 



X 



< 1 



and 



X* 



Poo 



o(x). 



(24) 



Second, the pressure at the point Q of the disturbed flow 
field should not be much smaller than the pressure imme- 
diately behind the shock in the intersection point of the 
shock surface with its normal through the point Q; i.e., 



= 0(1) (on x = const, y > 0) . 



(25) 



We follow Schneider (1968) to express the pressure in an 



(12) arbitrary point Q approximately by (see Appendix) 



P = P 



K + 2 [^* - KP, S*)]} 1/2 oV • (26) 



It may be noted that all quantities on the equation's right- 
hand side are given by the boundary conditions at the shock 
or by the equation of state. Terms coming from the excluded 
portions of the flow field - namely the stagnation region, as 
well as the region near the stagnation region where u <^u- 
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are of the order of x and therefore contribute only negligibly 
to the integral in equation (26). Consequently, the whole 



gas-dynamic state is known in the streamline coordinate 
system (x, ip) by evaluating S = S*(ip) and P from equation 
(26). 

The actual location of the body in space can be deter- 



mined by solving the differential equation (18) via separa- 



tion of variables giving the distance from the shock surface 
y as a function of x and 

l_ /^oV 



y 



j cos j3 



2r 



-y 



(27) 



Neglecting errors of 0(x) we may replace the velocity com- 
ponent in x-direction by 

u 2 = ul + 2 [ft* - ft(P, S*)] + . . • , (28) 



which follows from the energy equation (13). The integral 
in equation (27) is then 



Y 



f 

Jib 



dV>' 



(29) 



H p(P,5*)K + 2 [ft*-ft(P,5*)]} 1 /2 

Solving the quadratic equation in y on the left-hand side of 
equation (27), we have to distinguish between plane (j = 0) 
and axisymmetric (j = 1) flows. Thus we have 

0: y = Y; 

2Y cos p\ 



for j 
for j 



1: y 



cos /3 



1 



(30) 
(31) 



To finally give these results in the convenient coordinates 
x and 2/, the transformations and Q have to be carried 
out. 

In the astrophysical context, the case of a perfect gas 
with constant specific heats is of great interes t. The in- 



verse co mpression ratio is then given by (see e.g. Landau & 
Lifshitzp987 ) 

= Til i 2 (w\ 

X * 7 + 1 (7 + l)M^sin 2 A/ [6Z> 

We point out that an analogous relation is valid for x, if all 
asterisks in equation (32) are replaced by hats. 

Using the shock conditions (fl9] )— ( 22 ), together with 
equation (32), the two integrals (|26p ana (29), which have 



to be evaluated, become 



Uook 



cos 2 p. 




Y = 



I 



cos /?* 



X*(Psin 2 p*/P sin 2 $) 



l/ 7 



7 _1 - 



Psin 2 $ 
Psin 2 /3* 



- sin 2 /?* 

-1/2 



(34) 



The curvature of a curve in space can be calculated by the 
well-known formula 



d z r 



dz 2 



3/2 



(35) 



In addition, the following simple relations for the the shock 
inclination angle turn out to be quite useful: 



dr 

tan 6 = — =: q , 
dz 

sin 2 p = 



cos 2 P 



l + q 2 ' 
1 



(36) 
(37) 
(38) 



1 + q 2 ' 

The density for a perfect gas is represented by the quantity 

-l 



P = P* 



P\ lh 



* \ P 



7 + 1 (7 + l)M^sin 2 & 



P sin 2 /3 
P sin 2 P* 



l/ 7 



(39) 



For simplicity we introduce dimensionless units; e.g. the 
normalized stream function then is \I/ = ip / PoqUoqL^ 1 , 
where L is a characteristic length. 

On the surface of the body we need ip = 0, so the pres- 
sure on the body surface P\>(x), as well as the shock layer 
thickness A(x), can be obtained by replacing the lower lim- 
its in equations ([33]) and ([34]) by zero: 



P h 



+ (1-X) sin 2 /3 




+ sin P* 



(40) 



where the relation 
P 



Poo poo U^. 

1 



has been used. Furthermore, we get 



- + (1 - X) sin 2 P 
+ (1 - x) sin 2 P 



(41) 



for j = : A = Y 



for j = l: A 



(42) 



<0=O 



cos P 



2r|^ =0 cos/3 



1/2' 



(43) 
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Fig. 2. Comparison between the numerical and analytical 
solution of a bow shock wave generated by a hypersonic 
flow (Mqo = 500) hitting the solid unit sphere. A per- 
fect monoatomic gas is assumed (7 = 5/3). The analyti- 
cal body prediction is represented by the thin solid line. 
The bow shock wave (thick solid line) is parametrized by 
r = (2.620 • £) " 471 with j = 1. The (dimensionless) density 
profile shown is colour-coded, with red the high and blue 
the low density. 

with 



Y 



X*(Psin 2 /3*/Psin^/3) 



2 mV7 



x < cos f3* + 



(7 - i)m; 



sin /?* 



(44) 



P sin 2 f3 
Psin 2 & 



-1/2 



d*. 



As soon as the bow shock wave is parametrized, the gas- 
dynamical quantities can be computed by evaluating two 
integrals (namely equations ( 40 ) and (|42| or ( 43 ) and ([44]) , 
together with the boundary values (|14|) and ( 19)-(p2[)). 



3. Results 

3.1. Validation 

The applicability of the analytical model for an astrophys- 
ical context is tested by simulating the bow shock wave 
induced by the interaction of an uniform hypersonic flow 



with the solid unit sphere. For the sake of simplicity we 
have chosen = Uoo = 500, which is ensured by setting 
Poo =7 = 5/3 and Poo = 1. A simulation snapshot over- 
laid by the analytical solution is depicted in Fig. [2] At the 
stage shown here, the simulation was evolved long enough 
(t = 0.2) so that the bow shock wave has already become 
stationary for some time. The modulus of the percentage 
error of the stagnation distance prediction is ~ 6.2 %, and 
the modulus of the maximum percentage error of the ob- 
stacle estimation is ~ 11.6%. These errors may arise from 
two sources; first, the analytical method is based on the 
assumption that the density ratio x across the shock is 
very small. Since we consider here a perfect monoatomic 
gas with 7 = 5/3, an error that can reach the 25% is in- 
troduced (cf. Eq. ([32])) . However, for the astrophysically 
more relevant case 01a cooling strong shock, the compres- 
sion ratio of four only represents a lower limit. Thus, when 
the cooling time proceeds, relation (24) becomes progres- 



sively more precise. Second, the hypersonic gas flow around 
a sphere features a complex wake that contains inter alia a 
second shock wave as a result of overexpansion with follow- 
ing compression. This explains especially the increasing de- 
viation of the analytical solution from the numerical result 
at (or in the vicinity of) the sphere's equator (Schneider, 
priv. comm.). 

3.2. Application to observations 

To model the scenario in SQ, where a large-scale bow 
shock wave is believed to be formed by the gas-rich spi- 
ral g alaxy NGC 7318b as i t plunges through 100 K HI 
gas ( |Trinchieri et aL 2003), the upstream Mach num- 
ber has chosen to be = 930. We assumed a plane 
(j = 0) supersonic flow of perfect monoatomic gas (7 = 
5/3) and parametrized the bow shock time-dependently 
(e.g. by f 2 = 1.81 z for t = 130 Myr after the start of the 
interaction). Results are plotted in Figs. [3] and [4j 

The standoff distance of the shock from the nose of the 
projectile (in this case from the leading edge of the galaxy 
NGC 7318b of SQ) can be immediately determined to be 
A ~ 2.3 x 10 22 cm (~ 7.4kpc). In Fig. analytically de- 
rived maps of several fluid variables for the galaxy-IGM 
interaction process in SQ are shown, namely pressure P 
(upper left panel), density p (upper right panel), tempera- 
ture T (lower left panel), and the approximate Mach num- 
ber (since we have ignored v in the derivation, which is 
the velocity component in ^-direction) u/c s ~ M (lower 
right panel), where c s is the local downstream velocity of 
sound. Since values of z that are too low lead to a singular- 
ity in the (numerical) integration, the plane and near plane 
post-shock region are missing. We have written a Fortran 
program based on the presented method that allows the 
hypersonic blunt body problem to be solved for the flow 
quantities in the whole post-shock flow field. Required in- 
put variables are the parametrized bow shock f (£), the flow 
geometry factor j, the upstream Mach number Moo, and the 
adiabatic index of the gas 7. Furthermore, upstream den- 
sity, length scale, and the velocity of the obstacle (or flow, 
respectively) have to be known to allow an easy comparison 
to observational data. This tool is available onlinqi] 

To test our analytical results we performed several two- 
dimensional numerical hydrodynamical simulations. We 
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Fig. 3. Upper panel Shape of the SQ galaxy NGC 7318b 
(thin line) and the bow shock (thick line) at t = 130 Myr 
(j = 0, f 2 = 1.81 z). A perfect gas with 7 = 5/3 and 
an upstream Mach number of = 930 have been as- 
sumed. Lower panel. Normalized pressure P/pooU^ along 
the galaxy surface (i.e. along the body streamline) (dashed 
line) and directly behind the bow shock wave (solid line). 



used the VH-1 hydrocode as a basis, which was written by 
the Numerical Astrophysics Gr oup at the Virginia Ins titute 
for Theoretical Astrophysics ( |Blondin et al. 1990). The 
code is based on a third-order accurate extension of the 
grid-based Godunov scheme, namely a Lagrangian-remap 
version of the piecewise parabolic method (PPMLR), and 
features thus good shock capturing. Optically thin radia- 
tive cooling via operator splitting has been included. Each 
numerical simulation has been carried out on an uniform 
2400 x 1200 Cartesian grid, with a computational domain of 
4.8 x 10 23 cm by 2.4 x 10 23 cm (~ 156 x 78kpc). This results 
in a grid resolution of Ax = Ay = 2 x 10 20 cm (~ 65 pc). 
Runs with other grid sizes showed that the output does 



not depend on resolution. The only difference is that small- 
scale instabilities and fragmentations of stripped material 
cannot be resolved on coarser grids. "Outflow" boundary 
conditions have been used along the top, the bottom, and 
on the right of the grid, while the left boundary is "in- 
flow" . The galaxy itself was modelled by a simple stratified 
gaseous ellipse. Each computation begins with the galaxy 
placed face-on at rest in an uniform supersonic stream of in- 
tergalactic gas, corresponding to the situation of the galaxy 
rushing through the IGM at constant speed. 

Figure [5] gives an impression of our simulations. In this 
example, we tried to reconstruct the conditions that pre- 
sumably occur in SQ. The ISM of the intruding galaxy 
NGC 7318b (whose semi-axes have been taken from the 
neeQ consists in our simple model of three equidistant 
layers, which are in pressure equilibrium; namely, an Hi 

3 , T core = 100 K), a Lockman 



core (/?core = 6.2 x 10 g cm 



layer (p Loc k 



-24 

1.0 x 10~ 25 



gem" 3 , T Lock = 6 000K), and a 



Reynolds layer (p Rey = 6.9 x 10~ 26 gem" 3 , T Rey = 8 000K). 
The IGM (pigm = 6.7 x 10- 27 gcm ~ 3 , T lGM = 100 K) 
flows from left to right at a constant hypersonic speed 
(Moo = 930). Solar metallicity is assumed in all the consid- 
ered media. In Figure [5] the colour-coded log temperature 
map (in K) with the superimposed velocity field is visual- 
ized at the times 60, 130, and 200 Myr after the start of 
the simulation. Owing to the reasonably good numerical 
resolution of the simulation the presented snapshots allow 
tracking of a multitude of hydrodynamical effects (such as a 
prominent bow shock wave, Kelvin-Helmholtz instabilities, 
Ray leigh- Taylor instabilities, gas compression in the course 
of radiative cooling, a turbulent Karman wake, etc.) that 
lead to a subsequent evolution of the galaxy's ISM. Please 
note that more sophisticated raagnetohydrodynamic simu- 
lations with an underlying galaxy dark matter halo (which 
are beyond the scope of this paper) would be required to 
study a possible reduction or even suppression of fluid dy- 
namical instabilities. However, at high Mach numbers, as 
they occur in the considered galaxy-IGM interactions, the 
velocity potential (kinetic energy density) V ~ u 2 /2 (with 
u ~ 1400 km/s) is at least five times higher than the grav- 
itational potential (including dark matter) $ ~ GM/R m [ n 
(with M ~ 2 x 10 45 g, and a minimum distance from the 
galaxy's centre of R m i n ~ 7x 10 22 cm). Thus, the bow shock, 
as well as the shock layer, are unlikely to be influenced sig- 
nificantly by gravity in general and dark matter in partic- 
ular. 

The calculated post-shock temperature at the intersec- 
tion between the galaxy's axis of symmetry and the tip of 
the bow shock is about 2.7 x 10 7 K (k B T ~ 2.33 keV). In 
contrast, at the bow shock wings, where the shock incli- 
nation angle reaches approximately 30°, the downstream 
temperature has dropped to about 6.3 x 10 6 K (/cb T ~ 
0.54 keV). Despite the simple assumptions, this value is 
in re markably good agreeme nt with the observational data 
since |Trinchieri et al. (2003) measured in SQ a post-shock 
gas temperature of 0.5 keV. 



3.3. Comparison of the analytical with the numerical model 

In Figs. [6] and [7] our analytical results are superimposed 
on the numerically derived maps. In the lower panels the 
percentage errors are given. We used a shock tracer that 
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Fig. 4. Analytically derived colour-coded maps of several flow quantities, namely log pressure P in dyncm -2 (upper left 
panel), log density p in gem -3 (upper right panel), log temperature T in K (lower left panel), and the velocity ratio 
u/c s ~ M (lower right panel) for the SQ galaxy NGC 7318b at the time t = 130 Myr (j = 0, f 2 = 1.81 z). A perfect gas 
with 7 = 5/3 and an upstream Mach number of = 930 have been assumed. 



scans the computational domain for entropy jumps in or- 
der to parametrize the numerically calculated bow shocks. 
In addition to numerical models of the galaxy-IGM inter- 
action in SQ, we also studied the galaxy group around 
IC 1262. Chandra and XMM-Newton observations have 
uncovered an extended cool X-ra y ridge at the centre of 
this group (Trinchieri et al. 2007). The authors speculate 



that this peculiar structure may nave been created by ef- 
ficient ram pressure stripping of the bright spiral member 
IC 1263 as it had passed through the group's dense core 



Oigm = 8.0 x 10- 27 gcm- 3 , Tigm = 2.0 x 10 7 K) at a rela- 
tive velocity of ~ 1000 km s -1 (translating to a Mach num- 
ber of about 5). Then, gas supersonically streaming behind 
the bow shock wave may have heated up the stripped ISM 
to the observed X-ray temperatures. A result of our numer- 
ical and analytical computation of this interaction is shown 
in the lower right panel of Fig. [6j 

It is no surprise that the errors grow with increasing 
distance from the stagnation point and time, since instabil- 
ities, which cannot be predicted analytically, begin to grow 
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Fig. 5. Colour-coded log temperature map (in K) overlaid 
with the velocity field of the SQ simulation at the times 
(a) t = 60Myr, (b) t = 130 Myr, and (c) t = 200 Myr. 
The vectors drawn show the direction and magnitude of 
the velocity field. 



owing to the increasing influence of shear flows. Moreover, 
it must be noted that mainly the dense and cold ISM phases 
of the galaxies (blue to light green regions in Figs. [6] and 
[7]) drive the bow shock waves. Therefore these regions are 
the ones enclosed by the analytical solution. Bearing these 
aspects in mind, the agreement is quite good for all con- 
sidered scenarios (even for the lower Mach number case). 
Furthermore, the analytically derived fluid variable maps 
(Fig. are in good agreement with the numerical calcula- 
tions (compare the shock layer temperature gradient shown 
in the lower left panel of Fig. [4] with the upper right panel 
of Fig. [6| . The modulus of the percentage error of the log 
temperature immediately behind and in the vicinity of the 
bow shock wave lies below 1 %. Also the detachment dis- 
tances are determined quite satisfactorily and the modulus 
of the percentage error ranges from ~ 5.9 % to ~ 38.5 %. A 
more detailed discussion of the flow and comparison to new 
detailed X-ray data of SQ and the group around IC 1262 is 
the subject of a forthcoming paper. 

Finally, it should be noted that the astrophysical pos- 
sibility of this analytical method is hardly restricted to 
galaxy-IGM interaction alone; on the contrary, any sce- 
nario, such as bow shocks around subclusters, gas bullets, 
jets (Herbig-Haro objects) or individual stars can be cov- 
ered with it as well. 



4. Conclusion 

This paper presented an analytical solution technique for 
the inviscid hypersonic blunt body problem and applied 
it for the first time to the interaction between galaxies 
and the IGM. The method follows an inverse approach, 
which means that the shape of the bow shock has to be 
parametrized to be consistent with the galaxy shape. Such 
an inverse approach has no substantial disadvantage in the 
view of the fact that at high impact velocities the shape of 
the bow shock becomes similar to the shape of the body. 
For lower Mach numbers we find that the problem can be 
solved iteratively. First, we assumed the geometric form of 
the bow shock, then analytically derived the shape of the 
galaxy, and compared it to the relevant parameters (e.g. 
axes ratio, dimensions). In the next step, the bow shock 
curve was changed, and the body shape recalculated until 
convergence was obtained. In all cases we tested, this was 
achieved fairly rapidly. The solution, which is valid in the 
whole flow field, and in particular takes velocity gradients 
into account along streamlines, is based on two main as- 
sumptions; first, the density ratio across the shock has to 
be large, and, second, the pressure at a point Q (Fig. [I]) of 
the disturbed flow field must not to be very small in com- 
parison to the pressure immediately behind the shock in 
the intersection point of the shock surface with its normal 
through Q. In our derivation, heat conduction, viscosity, as 
well as terms of the order of \ (density ratio) are neglected. 
It should be stressed that it is not required that the shock 
layer, which is the area between the bow shock wave and the 
projectile, has to be thin. The presented treatment thus sur- 
passes in sophistication methods like the analytic solution 
for the thin-shell probl em of a hypersonic w ind interacting 



with a rigid sphere by |Canto fc Raga (1998) and allows as- 
trophysical problems that involve bow shocks generated by 
a hypersonic flow to be tackled more realistically. 

For a symmetrically assumed galaxy at a zero angle 
of attack, the stagnation streamline and the stagnation 
point are along the centreline. The stagnation streamline 
crosses the bow shock at precisely the normal shock point 
(/? = 7r/2) and hence the entropy of the stagnation stream - 
line in the shock layer is maximal (Hayes & Probstein 1966 ). 
However, we found that the galaxy need not necessarily be 
axially symmetric oriented with respect to the flow direc- 
tion. Then the stagnation streamline does not pass through 
the normal portion of the shock wave so it does not co- 
incide with the maximum entropy streamline that always 
intersects the shock at right angles. For all possible (two- 
dimensional) galaxy inclination angles, however, these two 
streamlines lie quite close to each other. Consequently, the 
presented method is applicable if the z-axis of the co- 
ordinate system is parallel to the velocity vector of the 
galaxy and if it goes through the stagnation streamline. 
The galaxy's bow shock halves above and below the z-axis 
then have to be evaluated separately. 

Our analytical investigation of galaxy-IGM interactions 
are in good agreement with results of the two-dimensional 
numerical simulations carried out with the (pure) hydro- 
dynamics code VH-1, which is PPMLR-based and therefore 
ideally suited to analysing shocks. As long as the initially 
made assumptions are not violated by Mach numbers that 
are too low, accurate predictions of the fluid quantities in 
the post-shock flow field, as well as the shock's stand-off 
distance, are possible. The high potential of the present 
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Fig. 6. Comparison of numerically (here log temperature in K is colour-coded) and analytically derived results. The thin 
line represents the outer envelope of the galaxy, whereas the thick line represents the bow shock. Upper left. Model of 
the interaction between NGC 7318b and the IGM in SQ at the time t = 60Myr (j = 0, f 2 = 1.82 z). Upper right. Same 
as upper left, but at at t = 130 Myr (j = 0, f 2 = 1.81 z). The lower panels give the corresponding percentage errors. 



analytical approximation suggests an application to many 
other astrophysical problems that involve bow shocks. 
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Appendix A: Derivation of the pressure in the 
shock layer 



The derivation of equation (26) in the main text is pre- 



sented here in detail. The momentum equation (15) can be 



10 



M.M. Schulreich & D. Breitschwerdt: An analytical solution for the hypersonic blunt body problem in the IGM 

AS 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 
z/ 10 23 cm 



Fig. 7. Like Fig.^ Upper left Model of the interaction between NGC 7318b and the IGM in SQ at the time t = 200 Myr 
(j = 0, f 21 = 1.40 z). Upper right. Model of the interaction between an IC 1263-like galaxy and the IGM in the galaxy 



group around IC 1262 at the time t = 130 Myr (j = 0, f 21 - 
The lower panels give the corresponding percentage errors. 



1.39 z). The upstream Mach number here is only = 5. 



rewritten as 



ld(u 2 
2 



dx 



Pd\nP 

p dx 



(A.l) 



From equations (19) and (20) we can estimate the ratio of 



the velocity components 

4=0(x 2 tan 2 /3) 



O 



(A.2) 



where and v\> are the velocity components at the surface 
of the body. Because the entropy increases with increasing 
shock inclination angle /3, it is valid to write S = (D(S*). 
Applying approximation (25), we have h = O(h) for the 



enthalpy. Moreover, it turns out that 
u = 0{u) , 



(A.3) 



where the energy equation (n3| has been used. In contrast, 
since the velocity increases along a streamline with decreas- 
ing pressure, it holds that 1/u = (D(l/u*). On the other 
hand, we can infer from the Rankine-Hugoniot jump condi- 
tions (19)-(22) that \ju — 0{\ju). This is true everywhere 



in the shock layer except the region, where cos f3* is very 
small compared to cos /?, i.e. where < ip < i/j. To include 
this region we have to write 



l-o ! , 

u \u h + [u-u h ]tp/tp 



(A.4) 



An analogous relation can be found for the density, if h and 
P in 



dp 



dp 
d~P 



dP 

h d ^ 



dp 

dh 



dh 

P dip 



(A.5) 



are replaced by means of equations (13) and fll6| ). As long 
as the thermodynamic functions h andP retain their order 
of magnitude, the thermodynamic functions (dp/dP)h and 
(dp/dh)p for a gas also do not change their orders of mag- 
nitude. The formal integration of the momentum equation 
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( 16 ) with the boundary condition P = P at ip = ip yields 

cos P 



F(P-P) 



l + O 



K + J- 



r 

* dv 
dx 



1 

lp — 1p Jrj) 



ydip 



(A.6) 



Furthermore keeping in mind that v is of the order of u 
everywhe re ap art from the st agnat ion region, we see from 
equation (A.5), together with (A.6), that the order of mag- 
nitude of dp/ dip does not depend on the stream function 
ip, except maybe in the region near the projectile surface, 
where u can be much smaller than u. Consequently, we ob- 
tain 

-=0\ 1 ,) . (A.7) 

P \p h + [p-p h }^/^J 



With the aid of our results ( |A.4 ), ( |A.7| ), and (11), equation 
(pb yields 



y = o 



In 



u h /u + [1 - u h /u)<p/<p \ 



v cos P u h /u - Pb/p^ p h /p+ [1 - p h /p\p/^p J ' 



Finally, combining equations (flTJ) and (|A.8|) gives 



O xsin/31n- 



(A.9) 



which can be combined with equations (21) and (25) to 
become 



^b 
^b 



Ph 



.1/2 



In 



Pb 



(A.10) 



The entropy remains constant on streamlines in the post- 
shock region. It is thus useful to introduce the effective 
isentropic exponent 



7e 



<91nP N 
dlnpy s 



lfdP\ 
P\dp)s 



(A.ll) 



which is not smaller than 1 for any gas. Consequently, the 
change in relative pressure along a streamline is not smaller 
than the change of relative density, i.e., 



dp < : 



O 



P 



p* 



If we now temporarily exclude the stagnation region (de- 
fined by t an 2 y ^> 1) from our consideration, equations 



( A. 2 ) and (A. 15) suggest neglec ting v 2 in comparison with 
u z in the momentum equation (A.l). Taking the isentropy 



of th e flow on streamlines into account, integration of equa- 
tion (A.l) yields 



InP 



InP* 



d(lnP) = 0. (A.16) 



s=s* 



If we once again exclude the stagnation region, as well as the 
region characterized by u <C ii, the first term of equation 
(A.16) can be rewritten by using 



0(u 2 ). 



(A.17) 



With the aid of equation ( A.13[ ) the integral in equation 
(A.16) can be evaluated: 



f ( 



d(lnP) = [&0 X *ln 



5=5, 



(A.18) 



Since the gas may expand substantially on the stagnation 
streamlines, the term ln(P*/P) can become very large. For 
further a nalysis it is advantageous to split the integral in 
equation (A.16) into two parts: 



An P 

! / - 

JlnP 



-) d(lnP) 

PJs=s* 



(A.19) 



d(lnP) = 0, 



s=s* 



where the second integral can be evalu ated b y using the 
assumption (25) together with equation (A. 13): 



AnP 
JlnP 



^)^d(lnP) = ^^ln^ 

= U 2 00 0( X *), 



(A.20) 



and thus turns out to be negligible. Equation (A.19) then 
simplifies to 



2 _ 2 
U (P) ~ U * 



dp 



s=s* 



0(X*U^ sin 2 P*) . (A.12) or, equivalently, 



On the other hand, the (dimensionless) quantity P/pooU^ 
cannot be very large for dp > 0, so for both dp < and 
dp > 0, we have 



\P) 



2[K -h{P,S*)\ 



(A.21) 



(A.22) 



P 

P 



(A.13) 



which also holds on the body streamline. Thus equation 
(A. 10) becomes 



where the well-known relation for the enthalpy dh = dP/ p 
(for dS = 0) has been used so that h = /i(P, S) is an 
equation of state of the gas. The subscript (P) implies 
that this approximation is used only t o calcu late the pres- 
sure P. Finally, by inserting equation (A.22), the Newton- 



v -± = o 



. Pb 

xx* — 
p , 



1/2 



ln^ 



Busemann pressure law can be rewritten to yield the im- 
portant result 



(A.14) 



P = P 



which can be simplified by using the assumption (24) and 
that p h /p = 0(1) 

— = 0(x). (A.15) 



{ul + 2[K - ft(P, S*)]} 1 / W , (26) 



used in Section [2j 



